Search for: All records

The radial distribution of gas within galactic haloes is connected to the star formation rate and the nature of baryon-driven feedback processes. Using six variants of the hydrodynamic simulation Simba, we study the impact of different stellar/AGN feedback prescriptions on the gas density profiles of haloes in the total mass range ${10}^{11}{\mathrm{M}}_{\odot }<M_{200\mathrm{c}}<{10}^{14}{\mathrm{M}}_{\odot }$and redshift interval $0<z<4$. We find that the radial profiles are well represented by a power law and that, for a fixed total halo mass, the slope and amplitude of such power law are generally weakly dependent on redshift. Once AGN-driven jets are activated in the simulation, the gas density profile of haloes with declines more gently with radial distance. We argue that this distinctive feature could be exploited with current observations to discriminate amongst the predictions of the different feedback models. We introduce a universal fitting formula for the slope and amplitude of the gas density profile as a function of total halo mass and redshift. The best-fit functions are suitable for all feedback variants considered, and their predictions are in excellent agreement with the numerical results. We provide the values of all fit parameters, making our fitting formula a versatile tool to mimic the effect of Simba feedback models onto N-body simulations and semi-analytical models of galaxy formation. Our results can also aid observational estimates of the gas mass within haloes that assume a specific slope for the underlying gas density profile. 

The distribution of offsets between the brightest cluster galaxies of galaxy clusters and the centroid of their dark matter distributions is a promising probe of the underlying dark matter physics. In particular, since this distribution is sensitive to the shape of the potential in galaxy cluster cores, it constitutes a test of dark matter self-interaction on the largest mass scales in the universe. We examine these offsets in three suites of modern cosmological simulations; IllustrisTNG, MillenniumTNG and BAHAMAS. For clusters above , we examine the dependence of the offset distribution on gravitational softening length, the method used to identify centroids, redshift, mass, baryonic physics, and establish the stability of our results with respect to various nuisance parameter choices. We find that offsets are overwhelmingly measured to be smaller than the minimum converged length scale in each simulation, with a median offset of in the highest resolution simulation considered, TNG300-1, which uses a gravitational softening length of . We also find that centroids identified via source extraction on smoothed dark matter and stellar particle data are consistent with the potential minimum, but that observationally relevant methods sensitive to cluster strong gravitational lensing scales, or those using the the “light traces mass” approach, in this context meaning gas is used as a tracer for the potential, can overestimate offsets by factors of $\sim 10$and $\sim 30$, respectively. This has the potential to reduce tensions with existing offset measurements which have served as evidence for a nonzero dark matter self-interaction cross section. 

ABSTRACT We investigate cosmological structure formation in fuzzy dark matter (FDM) with the attractive self-interaction (SI) with numerical simulations. Such a SI would arise if the FDM boson were an ultra-light axion, which has a strong CP symmetry-breaking scale (decay constant). Although weak, the attractive SI may be strong enough to counteract the quantum ‘pressure’ and alter structure formation. We find in our simulations that the SI can enhance small-scale structure formation, and soliton cores above a critical mass undergo a phase transition, transforming from dilute to dense solitons. 

ABSTRACT Bose–Einstein condensate dark matter (BECDM, also known as fuzzy dark matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established cold dark matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts (z ≳ 5) for a boson mass m = 2.5 × 10−22 eV, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered – a warm dark matter-like (‘WDM’) model often used as a proxy for BECDM. Our simulations confirm that ‘WDM’ is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/‘WDM’ compared to the CDM case: in BECDM/‘WDM’ first stars form at z ∼ 13/13.5, while in CDM star formation starts at z ∼ 35. The signature of BECDM interference, not present in ‘WDM’, is seen in the evolved dark matter power spectrum: although the small-scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM.